use "${filedata}for_regs.dta", clear
est clear
set scheme s2mono

/////DISTRIBUTION OF HARSHNESS figure
histogram iv_initial if tag_within_art==1 , ///
graphregion(color(white)) fcolor(emidblue)  ///
xtitle("Examiner Harshness") ytitle("Density")  
graph export ${results}exm_harshness.pdf, replace 

/////DISTRIBUTION OF HARSHNESS RESID figure
reghdfe iv_initial , abs(fe2=art_year) resid
predict IV2,resid
histogram IV2 if tag_within_art==1, ///
xtitle("Residualized Examiner Harshness") ytitle("Density")  w(0.04)
graph export ${results}exm_harshness_resid.pdf,replace 
